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During  this  period  our  program  for  evaluating  high-order  perturbation 
theory  was  essentially  completed,  and  significant  progress  was  made  in  direct 
numerical  solution  of  the  time-dependent  Schrodinger  equation.  In  addition, 
contributions  were  made  in:  understanding  the  role  of  spatial  dimensionality  in  the 
solution  of  model  systems  (e.g.  the  one-dimensional  delta-function  potential);  the 
development  of  a  general  R-matrix  code  for  perturbative  computations  of 
multiphoton  processes  in  many-electron  atoms;  and  multiphoton  ionization  with 
rvo  commensurate  laser  frequencies. 

The  main  development  in  the  work  on  high-order  perturbation  theory  invol¬ 
ved  the  extension  of  our  Sturmian  function  approach  to  treat  final  states  in  the  pho- 
toionization  continuum.  This  was  done  via  a  complex-coordinate-rotation  tech¬ 
nique,  which  gives  us  a  perturbative  expansion  for  the  complex  energies  of  a  non- 
Hermitean  Hamiltonian  in  powers  of  the  external  field  strength.  In  this  picture, 
all  atomic  states  become  resonances  in  the  presence  of  the  field;  the  imaginary  part 
of  the  complex  energy  of  a  resonance  is  inversely  proportional  to  its  lifetime.  The 
coefficients  in  the  expansion  are  independent  of  the  complex-coordinate  rotation 
angle,  which  is  a  significant  advantage  over  previous  methods.  Using  this 
approach,  we  are  able  to  compute  the  energy  shifts  associated  with  above- threshold 
ionization  (ATI)  channels,  in  which  photons  are  absorbed  by  a  continuum  electron. 
We  have  tound  that  the  onset  of  ATI  can  either  increase  or  decrease  the  total 
ionization  rate,  in  a  manner  which  depends  in  a  complicated  fashion  upon  the 
radiation  frequency.  This  is  differs  from  the  results  of  the  model  of  Deng  and 
Eberly,  in  which  ATI  always  acts  to  decrease  the  total  ionization  rate.  Comparisons 
made  with  independent,  non-perturbative  calculations  by  Chu  and  Cooper  showed 
that  corrections  to  lowest-order  perturbation  theory  can  be  significant  in  intensity 
regimes  where  perturbation  theory  is  effectively  convergent,  so  the  utility  of  this 
approach  is  not  limited  to  cases  in  which  lowest-order  perturbation  theory  is 
dominant.  However,  the  range  of  validity  of  perturbation  theory,  as  determined  by 
the  convergence  properties  inferred  from  numerical  calculations,  is  much  lower 
than  has  been  expected:  for  hydrogen  in  a  Nd:YAG  field  the  highest  intensity  at 
which  perturbation  theory  converges  is  of  the  order  3  x  1()D  W  cm~2.  The  nonlinear 
susceptibilities  that  describe  high-harmonic  generation  were  also  calculated  by  this 
method.  The  qualitative  behavior  of  the  high-order  susceptibilities  changes  when 
the  highest  virtual  state  lies  in  the  continuum  vs.  the  discrete  spectrum:  the 
critical  intensity  Ic,  at  which  perturbation  theory  breaks  down,  decreases  with 
increasing  order  when  ionization  is  energetically  forbidden,  but  exhibits  an  increase 
as  the  ionization  threshold  is  crossed.  These  results  have  been  published  in  the 
publications  numbered  4-6  below,  and  have  been  presented  in  invited  (1,2,4)  and 
contributed  (4)  talks  at  conferences. 

Perturbation  theory  was  also  applied  to  the  excitation  of  hydrogen  by 
radiation  at  two  frequencies,  to  and  2tu,  with  a  fixed  phase  relationship.  This  is  a 
problem  that  has  recently  become  of  experimental  interest.  Our  results  show  that 
the  ionization  rate  depends  upon  the  phase,  in  a  manner  that  is  determined  by  the 
conrir....,^  nhr.  ;e  shifts  of  the  electronic  wavefunction.  This  is  a  completely 
different  mechanism  than  the  enhanced  tunneling  that  is  predicted  in  classical 


treatments  of  the  problem.  This  work  was  reported  in  a  contributed  talk  (2),  but  has 
not  yet  been  submitted  for  publication. 

Work  continued  on  a  collaborative  project  with  K.  T.  Taylor  and  M.  Smith  at 
the  University  of  London,  in  which  a  general  R-matrix  code  is  being  developed  for 
the  computation  of  multiphoton  processes  in  many-electron  atoms  within  the 
framework  of  lowest-order  perturbation  theory.  This  code  is  nearly  complete,  and 
the  work  done  in  this  reporting  period  consisted  of  making  extensive  comparisons  of 
its  results  for  hydrogen  (e.g.  seven  photon  ionization)  with  those  obtained  by  the 
Sturmian  function  method.  Relatively  good  agreement  (at  the  10%  level)  has  been 
obtained,  except  near  intermediate  resonances,  whose  energies  are  slightly  different 
in  the  two  methods.  This  work  was  reported  in  a  contributed  talk  (3),  and  is  being 
written  up  by  M.  Smith  as  his  Ph.D.  thesis. 

The  frequency-dependent  polarizability  of  an  electron  bound  by  a  short  range 
potential  was  investigated  in  a  general  framework,  which  allowed  the  spatial 
dimensionality  of  the  system  to  be  treated  as  a  continuously  variable  parameter. 
This  was  motivated  by  the  frequent  use  of  delta-function  potentials  for  modelling 
the  response  of  atoms  to  strong  radiation  fields.  It  was  found  that  a  delta-function 
potential  taken  to  have  the  same  binding  energy  as  the  hydrogen  atom  has  a 
polarizability  much  lower  than  that  of  the  atom,  even  in  cases  far  from  resonance. 
The  polarizability  exhibits  a  monotonic  decrease  with  increasing  spatial 
dimensionality.  Thus  short-range  potentials  will  tend  to  underestimate  the 
response  of  a  real  atom  on  the  rising  edge  of  a  laser  pulse.  This  work  was  published 
(1)  and  presented  in  a  contributed  talk  (1). 

Numerical  integration  of  the  time-dependent  Schrodinger  equation  for 
hydrogen  in  a  radiation  field  began  in  earnest,  building  upon  preliminary  work  done 
late  in  1989.  The  main  advance  involved  the  adoption  of  massively  parallel 
processing  techniques  on  a  Connection  Machine  computer.  This  approach  divides 
configuration  space  into  many  (up  to  several  thousand)  distinct  regions,  each  of 
which  is  handled  by  a  dedicated  central  processing  unit  (CPU)  with  a  small 
amount  of  local  memory.  Each  CPU  integrates  the  equations  of  motion  within  its 
own  region,  using  the  contents  of  its  own  local  memory  and  that  of  its  nearest 
neighbors.  Thus  the  time  taken  to  solve  the  equations  is  roughly  independent  of  the 
number  of  regions,  if  a  sufficient  number  of  processors  are  available.  We  have 
integrated  the  equations  of  motion  for  hydrogen  in  radiation  fields  of  various 
frequencies  and  intensities  (up  to  1015  W  cm'2),  over  about  ten  optical  cycles. 
Machine  performance  of  several  gigaflops  has  been  attained,  which  is  still  somewhat 
below  the  theoretical  limit  for  the  configuration  we  are  using  (the  CM-2  facility  at 
the  Northeast  Center  foi  Parallel  Architectures,  Syracuse  University).  We  find 
the  plateau  phenomenon  in  the  harmonic  radiation  spectrum,  and  get  fieM-induced 
shifts  of  above-threshold  ionization  peaks  that  are  in  good  agreement  with  recent 
experimental  results  from  the  University  of  Bielefeld.  Work  has  been  started  on 
the  dynamics  of  two- electron  atoms  iu  strong  radiation  fields.  Our  presentation  of 
this  work,  which  is  still  in  a  developmental  phase,  is  limited  to  a  singte  invited  talk 
(7).  A  significant  amount  of  analytical  work  has  been  done  in  support  of  the 
numerical  effort.  This  has  led  to  a  comprehensive  understanding  of  the  propagation 
of  a  free  wavepacket  on  a  discrete  mesh,  and  to  a  closed-form  solution  of  the  bound  s 


states  of  hydrogen  in  the  discretized  representation.  This  has  been  presented  in  an 
invited  talk  (3),  and  is  being  written  up  for  publication. 

The  personnel  employed  on  this  project,  and  their  sources  of  support,  were  as 
follows: 


T.  J.  Mcllrath,  co-principal  investigator 
C.  W.  Clark,  co-principal  investigator 
L.  Pan,  research  associate 
J.  Parker,  research  associate 
S.  Blodgett-Ford,  graduate  student 


1  month  (AFOSR) 
3  months  (NIST) 

12  months  (AFOSR) 
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